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Abstract 

A generic feature of viable F{R) gravity is investigated: It is demonstrated that during the 
matter dominated era the large frequency oscillations of the effective dark energy may influence 
the behavior of higher derivatives of the Hubble parameter with the risk to produce some singular 
unphysical solutions at high redshift. This behavior is explicitly analyzed for realistic F{R) models, 
in particular, exponential gravity and a power form model. To stabilize such oscillations, we 
consider the additional modification of the models via a correction term which does not destroy the 
viability properties. A detailed analysis on the future evolution of the universe and the evolution 
history of the growth index of the matter density perturbations are performed. Furthermore, 
we explore two applications of exponential gravity to the inflationary scenario. We show how it is 
possible to obtain different numbers of e-folds during the early-time acceleration by making different 
choices of the model parameters in the presence of ultrarelativistic matter, which destabilizes 
inflation and eventually leads to the exit from the inflationary stage. We execute the numerical 
analysis of inflation in two viable exponential gravity models. It is proved that at the end of the 
inflation, the effective energy density and curvature of the universe decrease and thus a unified 
description between inflation and the ACDM-like dark energy dominated era can be realized. 

PACS numbers: 04.50.Kd, 95.36. +x, 98.80.-k 
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I. INTRODUCTION 



The current cosmic acceleration is supported by various observations such as Supernovae 
la (SNe la) jl|, large scale structure (LSS) ^ with baryon acoustic oscillations (BAO) jsl, 
cosmic microwave background (CMB) radiation 4-6| and weak lensing 3|- There exist 
two representative procedures to solve this problem, namely, introducing "dark energy" in 
general relativity (for recent reviews in terms of dark energy, see SMlOl]) and modifyingthe 

3)- 
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gravitational theory like F{R) gravity (for recent reviews on modified gravity, see 
In this paper, we adopt modified gravity approach to describe the inflation and dark energy 
eras. 

There proposed several viable F{R) gravity models have been constructed (for concrete 



viable models, see, e.g., the above reviews or 



16| and references therein). The conditions for 



the viability are summarized as follows: (i) Positive definiteness of the effective gravitational 
coupling, (ii) Matter stability condition 17H19|. (iii) In the large curvature regime, the 
model is close to the A-Cold-Dark-Matter (ACDM) model asymptotically, (iv) Stability 



of the late-time de Sitter point 

tests 



17 



22 



20 



21| . (v) The equivalence principle, (vi) Solar-system 



23| . It is considered to be one of the most significant issues on cosmology in the 
framework of F(,R) gravity to realize the unification of inflation with the late time cosmic 



acceleration 
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dark energy, see 



24 



17 1 (for the first proposal oi 1/R gravity as gravitational alternative for 



25|). 



In this paper, we study a generic feature of viable F{R) gravity models, in particular, 
exponential gravity and a power form model. We find that the behavior of higher derivatives 
of the Hubble parameter may be influenced by large frequency oscillations of effective dark 
energy, which makes solutions singular and unphysical at a high redshift. Therefore, in order 
to stabilize such oscillations, we examine an additional correction term to the model and 
remove such an instability with keeping the viability properties. We also demonstrate the 
cosmological evolutions of the universe and growth index of the matter density perturbations 
in detail. Furthermore, by applying two viable models of exponential gravity to inflationary 
cosmology and executing the numerical analysis of the inflation process, we illustrate that 
the exit from inflation can be realized. Concretely, we demonstrate that different numbers 
of e-folds during inflation can be obtained by taking different model parameters in the 
presence of ultrarelativistic matter, the existence of which makes inflation end and leads to 
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the exit from inflation. Indeed, we observe that at the end of the inflation, the effective 
energy density as well as the curvature of the universe decrease. Accordingly, a unified 
description between inflation and the late time cosmic acceleration is presented. We use 
units of /cb = c = ^ = 1 and denote the gravitational constant SttG by = 87r/Mpi^ with 
the Planck mass of Mpi = G'^/^ = 1.2 x lO^^GeV. 

The paper is organized as follows. In Sec. II, we briefly review the formulations of 
F{R) gravity. We use the fluid representation of F{R) gravity 26||. Here, in the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) background, the equations of motion with the addition 
of an effective gravitational fluid are presented. In Sec. Ill, we explain two well-known viable 
F{R) gravity models and show those generic feature occurring in the matter dominated era, 
when large frequency oscillation of dark energy appears and influences on the behavior 
of higher derivatives of the Hubble parameter in terms of time with the risk to produce 
some divergence and to render the solution unphysical. Thus, we suggest a way to stabilize 
such oscillations by introducing an additive modification to the models. We also perform a 
numerical analysis of the matter dominated era. In Sec. IV, we demonstrate that the term 
added to stabilize the dark energy oscillations in the matter dominated epoch does not cause 
any problem on the viability of the models, which satisfy the cosmological and local gravity 
constraints. We investigate their future evolution and show that the effective crossing of the 
phantom divide, which characterizes the de Sitter epoch, takes place in the very far future. 
We also analyze the growth index using three different ansatz choices. The second part of 
the paper is devoted to the study of F{R) models for the unification of the early-time cosmic 
acceleration, i.e., inflation, and the late-time one. In Sec. V, we explore two applications of 
exponential gravity for inflation. In particular, we show how it is possible to obtain different 
numbers of e-folds during inflation by making different choices of model parameters in the 
presence of ultrarelativistic matter in the early universe. In Sec. VI, we execute the numerical 
analysis of inflation and illustrate that at the end of it the effective energy density and the 
curvature decrease and eventually the cosmology in the ACDM model can follow. Finally, 
the summary and outlook for this work are given in Sec. VII. For reference, we also explain 
the procedure of conformal transformation in Appendix A and asymptotically phantom or 
quintessence modified gravity in Appendix B. 
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II. F{R) GRAVITY AND ITS DYNAMICS IN THE FLRW UNIVERSE: GEN- 
ERAL OVERVIEW 



In this section, we briefly review formulations in F{R) gravity and derive the gravitational 
field equations in the FLRW space-time. The action describing F{R) gravity is given by 

-F{R) 



M 



2k2 



^(matter) 



(2.1) 



where F{R) is a generic function of the Ricci scalar R only, g is the determinant of the metric 
tensor g^^^ ^(matter) jg ^.j^g matter Lagrangian and M. denotes the space-time manifold. In 
a large class of modified gravity models reproducing the standard cosmology in General 
Relativity (GR), i.e., F{R) = R, with a suitable correction to realize current acceleration 
and/or inflation, one represents 

F{R) = R + f{R). (2.2) 

Thus, the modification of gravity is encoded in the function f{R), which is added to the 
classical term R of the Einstein-Hilbert action in GR. In what follows, we discuss modified 
gravity in this form by explicitly separating the contribution of its modification from GR. 
The field equation simply reads 



^g,, {F{R) - RF'iR)) + (V^V. - g,,D) F\R) 



(2.3) 

Here, is the covariant derivative operator associated with g^y, = g'^'^'V u4> is the 
covariant d'Alembertian for a scalar field </>, and j'^Jf^**'''') = diag (p^, -Pm, -Pm) is the 
contribution to the stress energy-momentum tensor from all ordinary matters, with pm and 
Pm being the energy density and pressure of matter, respectively. Moreover, the prime 
denotes the derivative with respect to the curvature R. The trace of equation (12.31) in GR 
is trivial as P = — ft;22^(mattcr) ^j^j^ ji(matter) ^gj^g the trace of the stress energy momentum 
tensor of matter. On the other hand, for P(P) gravity, it becomes 



where 

dVeS _ 1 

dF'{R) 3 



(2P(P) - PP'(P) + K^T^"^^"'^'')) , (2.5) 
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with F'{R) being the so-called "scalaron", i.e., the effective scalar degree of freedom. On 
the critical points of the theory, the effective potential Ves has a maximum (or minimum), 
so that nF'{R) = 0. 

The most general flat FLRW space-time is described by the metric 

ds^ = -Af\t)dt^ + a{tfdi^ , (2.6) 

where a{t) is the scale factor of the universe and M{t) is an arbitrary function of the cosmic 
time t. In the following, we take the gauge Mit) = 1 and hence the Ricci scalar reads 

R = 12H^ + 6H, (2.7) 

where H = a{t)/a(t) is the Hubble parameter and the dot denotes the time derivative of 
dt{= d/dt). In the fiat FLRW background, from the (//, i/) = (0,0) component and the 
trace part of (/i, z/) = (z,j) (with z,j = 1, ■ ■ ■ ,3) components in Eq. (12. 3p . we obtain the 



gravitational field equations [11 1 



Peff = ' (2-8) 

P,s = -^(2H + 3H^) . (2.9) 



k2 . 

Here, p^s and Pgff are the effective energy density and pressure of the universe, respectively, 
defined as 

^ {F'R - F) - 6H^{F' - 1) - 6HF' , (2.10) 



PcflF = Pm + TT-;^ 
-PcflF = -Pm + 



{F'R-F) + {4H + 6H^){F' -1)+4HF' + 2F' . (2.11) 



In this way, we have a fluid representation of the so-called geometrical dark energy in F{R) 
gravity with the energy density pde = PcS — P and pressure Pde = -PcflF — -P- However, it is 
important for us to remember that gravitational terms enter in both left and right sides of 
Eqs. (12. 8p and (12. 9p . For general relativity in which F{R) = R, pcs = Pm and Pcfi = Pm and 
therefore Eqs. (12. 8p and (12. 9p lead to Friedman equations. 

We also explain basic equations that we use to carry out our analysis. In order to study 
the dynamics of F{R) gravity models in the fiat FLRW universe, we may introduce the 



variable [16i |27| 



yniz) ^P^ = ^-{z + lf-x{z + 1)' • (2.12) 
Pm(o) 



Here, Pm(o) is the energy density of matter at the present time, is the mass scale, given 
by 



and X is defined as 



a 



3 

Pr(0) 
Pm(0) 



3.1 X 10" 



where pi.(o) is the current energy density of radiation and 2 = l/a(t) — 1 is the redshift. Here, 
we have taken the current value of the scale factor as unity. By using Eqs. f l2.8p and fl2.12p . 
we find 

(fyH{z) , ^ dyniz) 



dz^ 



+ Ji- 



where 



Jl 



(z + l) 

1 



-3 



dz 



1 



J2{yH{z)) + J3 = 0: 



1 - F'{R) 



yH + {z + 1)3 + x{z + 1)^ Qfh?F"{R) 
1 2 - F'{R) 

yn + iz + 1)3 + x{z + 1)^ 3m^F"{R) 



[z + l)^ 
-3{z + l) 

(1 - F'{R)){{z + 1)3 + 2x{z + I f) + {R- F{R))/{3m^) 



iz + l)^{yH + {z + ir + x{z + ir) 
Furthermore, the Ricci scalar is expressed as 



Qm^F"{R) 



(2.13) 

(2.14) 
(2.15) 

(2.16) 



R = 3m 



Z.2 



4yH{z)-{z + l)^y^ + {z + lf 



(2.17) 



In deriving this equation, we have used the fact that —{z + l)H{z)d/dz = H{t)d/d{\na{t)) = 
d/dt, where H could be an explicit function of the red shift as H = H{z), or an explicit 
function of the time as H = H{t). In general, Eq. (12.131) can be solved in a numerical way, 
once we write the explicit form of an F{R) gravity model. 

III. GENERIC FEATURE OF REALISTIC F{R) GRAVITY MODELS IN THE 
MATTER DOMINATED ERA 



In this section, we consider viable F{R) gravity models representing a realistic scenario 



to account for dark energy, in particular, two well-known ones proposed in Refs. |27H30| (for 
more examples and detailed explanations on viable models, see, e.g., 16|, IMI] and references 
therein). We show that for these models, large frequency oscillation of dark energy in the 
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matter dominated era appears, and that it may influence on the behavior of higher derivatives 
of the Hubble parameter with respect to time. Such a oscillation has the risk to produce 
some divergence, and therefore we suggest a way to stabilize the frequency oscillation by 
performing the subsequent numerical analysis. In these models, a correction term to the 
Hilbert-Einstein action is added as F{R) = R+f{R) in fl2.2p . so that the current acceleration 
of the universe can be reproduced in a simple way. Namely, a vanishing (or fast decreasing) 
cosmological constant in the fiat limit of i? — )■ is incorporated, and a suitable, constant 
asymptotic behavior for large values of R is exhibited. 



A. Realistic F{R) gravity models 



Ref. 



32|) 



First, we explore the Hu-Sawicki model [27[ (for the related study of such a model, see 



C2(i?/m2)n + l ^ C2 C2(i?/m2)n + 1 ' ^-^-^^ 

where is the mass scale, ci and C2 are positive parameters, and n is a natural positive 
number. The model is very carefully constructed such that in the high curvature regime, 
m^Ci/c2 = 2A can play a role of the cosmological constant A and thus the ACDM model 
can be reproduced. 



Moreover, in Ref. 28| another simple model which may easily be generalized to reproduce 



also inflation has been constructed 

F{R) = R-2A[l- e-^/(''^)] , (3.2) 

where 6 > is a free parameter. Also in this model, in the flat space the solution of the 
Minkowski space-time is recovered, while at large curvatures the ACDM model is realized. 
This kind of models can satisfy the cosmological and local gravity constraints. Both of these 
models asymptotically approach the ACDM model in the high curvature regime. Indeed, 
however, the mechanisms work in two different manners, i.e., via a power function of R (the 
first one) and via an exponential function of it (the second one). For our treatment, we 
reparameterize the model (13.11) by describing cifn? /c2 = 2A and (02)^^" m? = bA with 6 > 0, 
so that we can obtain 



Through this procedure, in both of these models the term 6 A corresponds to the curvature 
for which the cosmological constant is "switched on" . This means 6 <^ 4, so that bA<^ 4A 
and hence R = 4A can be the curvature of de Sitter universe describing the current cosmic 
acceleration. In the mode in Eq. (13. 3p . since n has to be sufficiently large in order to 
reproduce the ACDM model, we have assumed n = 4 and we keep only the parameter b free. 



B. Dark energy oscillations in the matter dominated era 



Despite the fact that the models in Eqs. fl3.2p and 03.31] precisely resemble the ACDM 
model, there is a problem that in the matter dominated era the higher derivatives of the 
Hubble parameter diverge and thus this can make the solutions unphysical. This problem 



;hese models. In order to under- 
for dark energy perturbations 



originates from the stability conditions to be satisfied by 
stand this point, we start from the equations found in Ref. 
in the matter dominated era. The critical points associated with the matter dominated era 
for a generic F{R) model have been investigated in Ref. 2l|. The resultant conditions are 
given by 



F'{R) = , pr = 0, R = 3H^ 
F{R) = RF'{R) . 



(3.4) 



As a consequence, from Eq. (I2.10p . we have pcs — Pm- Suppose that yniz) ^ (1 + z)^ 
and x(l + z)^ <^ (-2 + 1)^ in order for dark energy and radiation to vanish in the matter 
dominated era, to the first order in terms of yniz), Eq. (I2.13P reads 

1 



yH{z)+y'Hiz) 



iz + l] 



[z + iy 



7 {1- F'{R))F"'{R) 
1 



RF"{R) 



+ 



2F"{Rf 

2(1 - F'{R))F"'{R) 



F"{Ry 



3 + 



2 - F'{R) - F{R)/R 



:z + i] 



(3.5) 



2RF"{R) 

where R = 3m{z + 1)^. In order to solve this equation locally, we can set z = zq + (z — zq), 
where \z — zo\ <^ zq, and perform a variation with respect to z. To the first order in terms 
of (z — zq), we find 



VHiz) + 



a P 

Vni^) + — —yniz) = Co + Ci{z - zo) 



{z - Zo) 



{z - ZoY 



(3.6) 



where 

7 (1 - F'{Ro))F'"{Ro) 



a 



2 2F"{Ro: 



2 



with and Ci being constants and -Rq = 3m^(zo + 1)'^. Thus, the solution of Eq. ( 13 .Gp is 
derived as 

yniz) = a + b ■ {z - zo) + Co ■ exp (^-a ± v^a^ _ (^ _ ^ (3_g) 

where a, b and Co are constants. Now, for the two models in Eqs. fl3.2p and (13. 3p . when 
i? > 6A, we find 

F'{R) ~ 1 , 

F"{R) ~ 0+ . (3.9) 

These behaviors guarantee the occurrence of the realistic matter dominated era. Further- 
more, since in the expanding universe {z — zq) < 0, it turns out that the dark energy 
perturbations in Eq. (13.80 remain small around Rq, and that we acquire 

{1- F'{Ro))F"'{Ro) _7 1 

2F"{Ror ^~2' RoF"{Ro) ^ ' ^^'^^^ 

for both these models. Owing to the fact that F"{R) is very close to 0^, the discriminant 

in the square root of Eq. ( 13.80 is negative and dark energy oscillates as 



3m^ 



yHiz) = ^+e 2(.o+i) 



A sin I V^^^'^ (_^ _ 2;g) ] + 5 cos I "^^^'"^ -(z — 20) 



(^0 + 1) J \izo + l) 

' (3.11) 

Here, A and B are constants and ai^2 and /3i^2 are given by Eq. (13. 7p . so they correspond 
to two models under investigation. In particular, ol\ = —3 for the model in Eq. (13. 2 p and 
a2 ~ -29/10 for the model in Eq. ((331), while /3i,2 ^ 1 / {RoF" (Rq)) , i.e.. 



/3i ~ — -— , (3.12) 



2i?o 

in case of exponential model in Eq. (13. 2 p and 



, M^Hm'i^y Ro fR^y 
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in case of model in Eq. (13. 3p . This means that the frequency of dark energy oscillations 
increases as the curvature (and redshift) becomes large. Moreover, the effects of such oscil- 
lations are amplified in the derivatives of the dark energy density, namely. 



OC (^"(^0))' 



(3.14) 



1 /2 

where J^{z) ~ (i? * F"{R))-''' /{z + 1) is the oscillation frequency and to is the cosmic time 
corresponding to the redshift zq. For example, we can analyze the behavior of the EoS 
parameter for dark energy defined as^ 

. ^ -Pde 



Pde 



3 yH[z) d[z) 



(3.15) 



By combining Eqs. (I3.1ip and (I3.15p . we find 

_ "1,2(^-^0) 

e 

0J-D^[z] ~ 



-1 + 



rn?e 2{zo+i) 



A 



A sin 



1,2 



(^0 + 1) 



[z — zq)\ + B' cos 



{zo + 1) 



(3.16) 



where 



A' 
B' 



■''-fA-^.B. 

"1,2 



(3.17) 



are constants. For large values of the redshift, the dark energy density oscillates with a 
high frequency and also its derivatives become large, showing a different feature of the dark 
energy EoS parameter in the models in Eqs. (13.20 and (13.30 compared with the case of 
the cosmological constant in GR. During the matter dominated era, the Hubble parameter 
behaves as 



{z + 1)=^/^ 



(3.18) 



2(;z + l)3/2_ 

If the frequency ^^(^o) in Eq. (I3.14p is extremely large, the derivatives of dark energy density 
could become dominant in some higher derivatives of the Hubble parameter as 



dV 



:H{t)\ 



to 



OC (Tizo))' 



> 1 



(3.19) 



which may approach an effective singularit y a nd therefore make the solution unphysical. 



We see it for specific cases. In Refs. 



16 



34 



35|, the cosmological evolutions in exponential 



^ Throughout this paper, we describe the EoS parameter by "w" and not "w" 
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gravity and the Hu-Sawicki model have carefully been explored. It has explicitly been 
demonstrated that the late-time cosmic acceleration which follows the matter dominated 
era can occur, according with astrophysical data. A reasonable choice is to take b = 1 for 
both these models. We also put A = 7.93m^ j^. We can solve Eq. f l2.13p numerically^ by 
taking the initial conditions aX z = Zi, where Zi ^ is the redshift at the initial time to 
execute the numerical calculation, as follows: 



d{z) 



0, 
A 



3m 



Here, we have used the fact that at a high redshift the universe should be very close to the 
ACDM model. We have set Zi = 2.80 for the model in Eq. (13.21) and Zi = 4.5 for the model 
in Eq. such that RF"{R) ~ 10"® at = 'ifn^{zi + 1)^. We note that it is hard to 

extrapolate the numerical results to the higher redshifts because of the large frequency of 
dark energy oscillations. 

Using Eq. (13.151) with yu, we derive cjde- In addition, by using Eq. (I2.17P we obtain R as 
a function of the redshift. We can also execute the extrapolation in terms of the behavior 
of Qde, given by 

fiDE(^) = ^ = ^ . (3.20) 

The numerical extrapolation to the present universe leads to the following results: For 
the model ([S2D, 2/h(0) = 2.736, wdeIO) = -0.950, ^deIO) = 0.732 and R{z = 0) = 
4.365, whereas for the model ([SJD, 2/h(0) = 2.652, wdeIO) = -0.989, 1^de(0) = 0.726 and 
R{z = 0) = 4.358. These resultant data are in accordance with the last and very accurate 

51, which are 



observations of our current universe 



'^de - -0.972+o;ogo , 

fioE = 0.721 ±0.015. (3.21) 



We have used Mathematica 7 ©. 
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Next, we introduce the deceleration g, jerk j and snap s parameters 36| 



^ ~ a{t) dt^ H{tf 

1 d^a{t) 1 _ H 

In what follows, we show the values of these cosmological parameters at the present time 
{z = 0) as the result of numerical extrapolation in our two models, which we called Model 
I in Eq. and Model II in Eq. (O, and the calculation in the ACDM model: 

q(z = 0) = -0.650 (ACDM) , -0.544 (Model I) , -0.577 (Model II) 

j{z = 0) = 1.000 (ACDM) , 0.792(ModelI) , 0.972 (Model II) 

s{z = 0) = -0.050 (ACDM) , -0.171(ModelI) , -0.152 (Model II) . 

The deviations of the parameters in Models I and II from those in the ACDM model are 
small at the present. However, since these parameters depend on the time derivatives of the 
Hubble parameter, it is interesting to analyze those behaviors at high curvature. Therefore, 
in Fig. [1] we plot the cosmological evolutions of q, j and s as functions of the redshift z. 
From this figure, we see that there exist overlapped regions for Models I and II with those 
in the ACDM model. 

The deceleration parameter in Models I and II remains very close to the value in the 
ACDM model, because in the first time derivative of the Hubble parameter the contribution 
of dark energy is still negligible. Hence, it guarantees the correct cosmological evolution of 
these models. However, it is clearly seen that in the jerk and snap parameters the derivatives 
of the dark energy density become relevant and the parameters grow up with an oscillatory 
behavior. Since the frequency of such oscillations strongly increases in the redshift, it is 
reasonable to expect that some divergence occurs in the past. We also remark that if from one 
side at high redshifts the exponential Model I is more similar to the ACDM model because 
of the faster decreasing of exponential function in comparison with the power function of 
Model II, from the other side it involves stronger oscillations in the matter dominated era. 

It may be stated that the closer the model is to the ACDM model (i.e., as much F"{R) is 
close to zero), the bigger the oscillation frequency of dark energy becomes. As a consequence. 
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(a) (b) 




(e) (f) 

FIG. 1: Cosmological evolutions of q{z) [(a) and (d)], j{z) [(b) and (e)] and s{z) [(c) and (f)] 
parameters as functions of the redshift z for Model I [(a)-(c)] and Model II [(d)-(f)] in the region 
of z > 0. 

despite the fact that the dynamics of the universe depends on the matter and the dark energy 
density remains very small, some divergences in the derivatives of the Hubble parameter can 
occur. In the models in Eqs. 03.21) and 03.31) . although the approaching manners to a model 
with the cosmological constant are different from each other, it may be interpreted that these 
models show a generic feature of realistic F{R) gravity models, in which the cosmological 
evolutions are similar to those in a model with the cosmological constant. The corrections 
to the Einstein's equations in the small curvature regime lead to undesired effects in the 
high curvature regime. Thus, we need to investigate additional modifications. 
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C. Proposal of a correction term 



In order to remove the divergences in the derivatives of the Hubble parameter, we in- 
troduce a function g{R) for which the oscillation frequency of the dark energy density in 
Eq. f l3.1ip acquires a constant value 1 / a/5, 5 > 0, for a generic curvature R ^ bA, and we 
stabilize the oscillations of dark energy during the matter dominated era with the use of a 

1 /3 

correction term. Since in the matter dominated era, i.e., z + 1 = [i?/(3m^)] , we have to 
require 

(3m^)^/3 ^ ^ 
R^/^ g"{R) ~ 6 

g{R) = -^A(—^j , 7>0, (3.23) 

where 7 = (9/2)5(3mV^) = 1-702 5. We explore the models in Eqs. ([321) and (JSJl) with 
adding these correction as 



R 



R 



F,{R) = R-2A{l-e'^)-^A(^] , (3.24) 



1/3 



F2{R) = R-2A ^ ~ ' ' 



1- ' 



We note that in both cases -^1,2(0) = and therefore we still have the solution of the flat 
space in the Minkowski space-time. The effects of the last term vanish in the de Sitter 
epoch, when R = 4L and these models resemble to a model with an effective cosmological 
constant, provided that 7 ^ (m^/A)-*^/^. We may also evaluate the dark energy density at 
high redshifts by deriving pde = Pcfi — Pm from Eq. f l2.10p and by putting R = 3m'^{z + 1)^ 
such that 

yH{z)^^[l + l{l + z)] . (3.26) 

According to the observational data of our universe, the current value of dark energy amount 
is estimated as yn = A/(3m) = 2.643. With the reasonable choice 7 ~ 1/1000, the effects 
of modification of gravity on the dark energy density begin to appear at a very high redshift 
(for example, at 2 = 9, yni^) = 1-01 x yniO)), and hence the universe seems to be very 
close to the ACDM model. However, while the pure models in Eqs. (13. 2p and (13. 3p mimic an 
effective cosmological constant, the models in Eqs. (I3.24p and (I3.25p mimic (for the matter 
solution) a quintessence fluid. Equation (I3.15P leads to 

+ 3(1 +(1+1)7) ' 
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so that when z ^ 7"^, wdeI-^) — —2/3. 

Thus, it is simple to verify that all the cosmological constraints 12[ are still satisfied. Since 
\F[ 2(-R ^ 6A) — 1| ^ 1, the effective gravitational coupling Gcfr = G iF[ ^{K) is positive, and 
hence the models are protected against the anti-gravity during the cosmological evolution 
until the de Sitter solution (-Rds = 4A) of the current universe is realized. Thus, thanks to 
the fact that |-F^"2(-R ^ &A) > 0|, we do not have any problem in terms of the existence 
of a stable matter. In Sec. llVt we also analyze the local constraints in detail, and we see 
that our modifications do not destroy the feasibility of the models in the solar system. It 
should be stressed that the energy density preserves its oscillation behavior in the matter 
dominated era, but that owing to the correction term reconstructed here, such oscillations 



keep a constant frequency T = a/1.702/7 and do not diverge. Despite the small value of 7, 
in this way the high redshift divergences and possible effective singularities are removed. 



D. Analysis of exponential and power- form models with correction terms in the 
matter dominated era 

In this subsection, we carry out the numerical analysis of the models in Eqs. fl3.24p and 
f l3.25p . In both cases, we assume 6 = 1 and 7 = 1/1000 and solve Eq. (I2.13P in a numerical 
way, by taking accurate initial conditions a.t z = Zi so that Zi ^ 2. By using Eq. (13. 26 p . we 
acquire 



dyniz) 



d{z) 



A 

3m 
A 

3m' 



where we have set Zi = 9. The feature of the models in Eqs. (13.240 and (I3.25P at the present 
time is very similar to those of the models in Eqs. (13.21) and (13. 3p . With the numerical 
extrapolation to the current universe, for the model in Eq. (I3.24p we have y_f/(0) = 2.739, 
ujbe{0) = -0.950, Qbe{0) = 0.732 and R{z = 0) = 4.369, while for the model in Eq. fl3:25|l . 
we find yniO) = 2.654, ujbe{0) = -0.989, ^deIO) = 0.726 and R{z = 0) = 4.361. We 
analyze those behaviors in the matter dominated era. It follows from the initial conditions 
yni^) = 2.670 and ci;de(9) = —0.997 that the universe is extremely close to the ACDM model 
also at high redshifts. We see how the dynamical correction of the Einstein's equation, which 
corresponds to, roughly speaking, the fact of having "a dynamical cosmological constant", 
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introduces an oscillatory behavior of dark energy density. Thanks to the contribution of 
the correction term, we obtain a constant frequency of such oscillations without changing 
the cosmological evolution described by the theory. In Fig. |21 we show the cosmological 
evolutions of the deceleration, jerk and snap parameters as functions of the redshift z in 
these models. There is overlapped region of the evolutions with those in the ACDM model. 
We may compare the graphics in Fig. [2] with the corresponding ones in Fig. [1] of the models 
in Eqs. (13. 2p and (13. 3 p without the correction term analyzed in Sec. Ill B. At high redshifts, 
the deceleration parameter is not influenced by dark energy and hence the behavior in both 
these models in Eqs. (I3.24p and (I3.25P are the same as that in the ACDM model. On the 
other hand, in terms of the jerk and snap parameters, the derivatives of dark energy density 
become relevant and accordingly these parameters oscillate with the same frequency as that 
of dark energy, showing a different behavior in comparison with the case of GR with the 
cosmological constant. However, here such oscillations have a constant frequency and do not 
diverge. The predicted value of the oscillation frequency is J-" = a/1. 702/7 = 41.255. The 
oscillation period is T = 271 / ~ 0.152. Thus, the numerical data are in good accordance 
with the predicted ones. (We can also appreciate the result by taking into account the fact 
that the number of crests per units of the redshift has to be 1/T ~ 7). 



IV. COSMOLOGICAL CONSTRAINTS AND FUTURE EVOLUTION 



In this section, first we show that the models in Eqs. (I3.24p and (13.251) satisfy the cos- 



mological and local gravity constraints 



23| . and that the term added to stabilize the dark 



energy oscillations in the matter dominated epoch does not cause any problem to these 
proprieties. The confrontation of F{R) models with SNIa, BAO, CMB radiation and grav- 
itational lensing has been executed in the past several works 37|. We have just seen that 
the models with the choice of 6 = 1 can be consistent with the observational data of the 
universe. Here, we examine the range of h in which the models are compatible with the 
observations and analyze the behavior of the models near to local (matter) sources in order 
to check possible Newton law corrections or matter instabilities. Then, we concentrate on 
the future evolution of the universe in the models and demonstrate that the effective crossing 
of the phantom divide which characterizes the de Sitter epoch takes place in the very far 
future. 
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(a) 



(b) 





(c) 




- 


(d) 




- 
- 





(e) (f) 
FIG. 2: Cosmological evolutions of q{z) [(a) and (d)], j{z) [(b) and (e)] and s{z) [(c) and (f)] 
parameters as functions of the redshift z for the model Fi(R) [(a)-(c)] and the model F2{R) [(d)-(f)] 
in the region of 2; > 0. 



In the way of trying to explain the several aspects that characterize our universe, there 
exists the problem of distinguishing different theories. It has been revealed that sometimes 
the study of the expansion history of the universe is not enough because different theories 
can achieve the same expansion history. Fortunately, theories with the same expansion 
history can have a different cosmic growth history. This fact makes the growth of the large 
scale structure in the universe an important tool in order to discriminate among the different 
theories proposed. Thus, the characterization of growth of the matter density perturbations 
become very significant. In order to execute it, the so-called growth index 7 



is useful. 



Therefore, in the second part of this section we study the evolution of the matter density 
perturbation for our F{R) gravity model. 
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A. Cosmological and local constraints 



We take 7 = 1/1000 in the models in Eqs. f l3.24p and f l3.25p . keeping the parameter b 
free. Now, the dark energy density is a function of z and b, i.e., yniz, b). We can again solve 
Eq. (I2.13P numerically, taking the initial conditions at = 9 as 



dyniz.b) 



d{z) 

yH{z,b) 



A 
A 



{l+^{Zi + l)) 



as we did in the previous section. We take 0.1 < 6 < 2. In Figs. [3] and HI we display the 
resultant values of dark energy EoS parameter oj-q^^z = 0, b) and flBE{z = 0, b) at the present 
time as functions of b for the two models. We also show the bounds of cosmological data in 
Eq. fl3.2ip . namely, the lines in rose denote the upper bounds, while the lines in yellow do 
the lower ones. By matching the comparison between the two graphics of every model, we 
find that in order to correctly reproduce the universe where we live with exponential gravity 
in Eq. (^M), 0.1 < b < 1.174, with power-law model in Eq. (K25^ . 0.1 < 6 < 1.699. The 
results are consistent with the choices in Sec. IIIIDI 





(a) (b) 

FIG. 3: Behaviors of cjdeC-z = 0,6) and of Qj)-e{z = 0,6) as functions of b for exponential model. 
The observational data bounds (horizontal lines) are also shown. 



Newton law corrections and stability on a planet surface 



In Ref. 39|, it has been shown that some realistic models of F{R) gravity may lead to 
significant Newton law corrections at large cosmological scales. We briefly review this result. 
From the trace equation (12. 4p . if we consider the constant background of i? = Rq, such that 
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(b) 



(a) 

FIG. 4: Behaviors of ix!de(-2 = 0,6) and of OdeC-z = 0,6) as functions of b for power-law model. 
Legend is the same as Fig. [3l 

2F{Rq) — RqF'{Rq) = 0, by performing a variation with respect to i? = i?^"^ + 5R and 
supposing the presence of a matter point source (like a planet), that is, T*^™^"™) = TqS{x), 
where 6{x) is the Dirac's distribution, we find, to first order in 6R, 

,2 



(□ - m^) 6R 



K 



To6{x) 



with 



The solution is given by 



m 



6R 



3F"(i?o) 

F'jRo) ^ 
3 \F"{Ro) ' 

-ToG{m^,\x\) 



K 



SF"{Ro 



where G{m?, \x\) is the correlation function which satisfies 

(□-m^) G(mMx|) = 5{x). 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



Hence, if m? < 0, there appears a tachyon and thus there could be some instability. Even if 
> 0, when is small compared with Rq, 6R 7^ at long ranges, which generates the 
large correction to the Newton law. For the pure exponential model in Eq. (13.21) without 
correction terms, when Rq ^ bA, m? reads 

(fe^A) «a 



2 

m ~ 



6 



(4.5) 



Therefore, in general m? / Rq is very large effectively. The same thing happens in the model 
in Eq. (13. 3p . Next, for the models in Eqs. (I3.24p and (I3.25P with correction terms, we have 



2A7 
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(4.6) 



Despite the fact that in this case is smaller than in Eq. (14. 5p . it still remains sufficiently 
large and the correction to the Newton law is very small. For example, the typical value of 
the curvature in the solar system is Rq ^ 10~^^eV^ (it corresponds to one hydrogen atom 
per cubic centimeter). In this case, from Eq. (14. 6 p we obtain m? / Rq ~ 2 x 10^. 

Concerning the matter instability [l^, this might also occur when the curvature is rather 
large, as on a planet {R ~ lO^^^eV^), as compared with the average curvature of the 
universe today [R ^ 10~^^eV^). In order to arrive at a stability condition, we can perturb 
again Eq. (12. 4p around R = Rb, where Rb is the curvature of the planet surface and the 
perturbation 6R is given by the curvature difference between the internal and the external 
solution. The curvature Rb = —ft;22^ (matter) (^gpg^cis on the radial coordinate r. By assuming 
6R depending on time only, we acquire 



- d1{5R) ^U{Rb)6R, (4.7) 



where 

U{Rb 



F"'{Rb)V F"'{Rb 



F"{Rb) J F"{R, 
F"'{Rb' 



g^^\^rRb^rRb ^ ' ^ 



3 SF"{Rb) 

{2F{Rb) - RbF\Rb) - Rb) . (4.8) 



^{F"{Rb)f 

Here, g^y is the diagonal metric describing the planet. If U{Rb) is negative, then the per- 
turbation 6R becomes exponentially large and the whole system becomes unstable. Thus, 
the planet stability condition is 

UiRb) > . (4.9) 

For our models in Eqs. (I3.24P and (13. 25 p . U{Rb) — m? , where m? is given by Eq. (14. 6 p 
again. Also in this case, we do not have any particular problem. For example, by putting 
Rb — 10~^*eV^, we find U{Rb)/Rb — 4 x 10^^. Thus, the models under consideration easily 
pass these local tests. 

B. Future universe evolution 

In de Sitter universe, we have R = R^s, where R^s is the constant curvature given by the 
constant dark energy density yn = Ho, such that yo = -Rds/12m^. Starting from Eq. (12.131) . 
we are able to study perturbations around the de Sitter solution in the models (I3.24p and 
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fl3.25p which provide this solution for i?ds = 4A and well satisfied the de Sitter condition 
2F(i?ds) = RdsF^'iRds) as a consequence of Eq. (12.41) in vacuum. Performing the variation 
with respect to yniz) = yo + yi{z) with |yi(2;)| ^ 1 and assuming the contributions of 
radiation and matter to be much smaller than yQ, at the first order in yi{z) Eq. (I2.13P reads 



+ 



a 



dyi{z) 



/3 



where 



a 



-2, /3 



dz {z + iy 

^F'iRds) 



y,iz)=Aaz + l), 
1 - F'iR^s) 



RF"{R, 



C = i + 



RdsF"{Rds] 



The solution of Eq. (I4.10p is given by 



yniz) = yo + yi{z) , 

y,{z) = Coiz + i)K-»±v/(^^^) + |(, + 1)3 , 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



where Co is a constant. The well-known stability condition for the de Sitter space-time, 
F'{Rds)/{{Rds)F"{Rds)) > 1, is also vahd. It has also been demonstrated that since in real- 



istic F 
25/16 



(R) gravity models for the de Sitter universe F"{R) 0+, F' (Rds) / {RdsF" (Rds)) > 



40| giving negative the discriminant of Eq. fl4.13p and an oscillatory behavior to the 



dark energy density during this phase. Thus, in this case the dark energy EoS parameter 
wde fl3.15p becomes 



<^de{R 



An cos 




(4.14) 



Bn sin 




RdsF"iRds: 



\og{z + 1 




and oscillates infinitely often around the line of the phantom divide wde = ~1 |4Q|. Ac- 
cording to various recent observational data, the crossing of the phantom divide occurred in 
the near past |4]|. These models possess one crossing in the recent past [16|, after the end 



of the matter dominated era, and infinite crossings in the future (for detailed investigations 



on the future crossing of the phantom divide, see 3l|), but the amplitude of such crossings 
decreases as (z + 1)^/^ and it does not cause any serious problem to the accuracy of the 



cosmologica! 



evolution during the de Sitter epoch which is in general the final attractor of 



the system 16|, |34]. However, the existence of a phantom phase can give some undesirable 



effects such as the possibility to have the Big Rip 



42| as an alternative scenario of the 
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universe (in such a case, the model may suddenly exit from ACDM description) or the dis- 
integration of bound structures which does not necessarily require to having the final (Big 



Rip) singularity 43|, . In this subsection, we show that in the models in Eqs. 



321 and 



fl3.25p the effective EoS parameter of the universe (for an alternative study, see 45|) defined 



as 



_Peff , , 2{z+l)dH{z) 
Wcff = ^ = -1 + oTj/ N -1 — 4.15 

never crosses the phantom divide line in the past, and that only when z is very close to 
— 1 (this means in the very far future), it coincides with ube and the crossings occur. We 
remark that Pcfi and Pcfr correspond to the total energy density and pressure of the universe, 
and hence that if dark energy strongly dominates over ordinary matter, we can consider 
WeflF ~ u^DE- In both of the models under investigation, we take again 7 = 1/1000 and keep 
the parameter b free, such that 0.1 < 6 < 1.174 (model in Eq. f l3.24p ) and 0.1 < 6 < 1.699 
(model in Eq. fl3.25p ). according to the realistic representation of current universe. The 
numerical evaluation of Eq. f l2.13p leads to H{z), given by 

H{z) = sjm^ \yi,{z) + (;2 + 1)3 + x{z + 1)^] , (4.16) 

and therefore ujc^^iz). We depict the cosmological evolution of as a function of the red 
shift z and the h parameter in Fig. [5] for the model in Eq. f l3.24p and in Fig. |6] for the model 
in Eq. f l3.25p . On the left panels, we plot the effective EoS parameter for — 1 < z < 2. 
We can see that for both of the models, independently on the choice of 6, cjde starts from 
zero in the matter dominated era and asymptotically approaches -1 without any appreciable 
deviation. Only when z is very close to —1 and the matter contribution to is effectively 
zero, we have the crossing of the phantom divide due to the oscillation behavior of dark 
energy. On the right panels, we display the behavior of the effective EoS parameter around 
z = —\. Here, we focused on the phantom divide line and we excluded the graphic area out 
of the range —1.0001 < Wcfr < —0.9999. The blue region indicates that Wcfr is still in the 
quintessence phase. We note that especially in the model in Eq. (13. 25 p . the first crossing of 
phantom divide is very far in the future. For example, with the scale factor a(t) = exp {HQt), 
where Hq ~ 6.3 x 10~'^^eV~^ is the Hubble parameter of the de Sitter universe, z = —0.90 
(when the crossing of the phantom divide may begin to appear in the exponential models) 
corresponds to 10^^ years. 
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(a) 



(b) 



FIG. 5: Cosmological evolution of lOcs as a function of the red shift z and the b parameter for the 
model in Eq. (j3.24p . The left panel plots it for — 1 < z < 2 and the right one displays around 



FIG. 6: Cosmological evolution of uj^s as a function of the red shift z and the b parameter for the 
model in Eq. p.25p . Legend is the same as Fig. [5l 

Avoidance of the phantom crossing with (inhomogeneous) fluid 

It may be of some interest to check if it is possible to avoid the crossing of the phantom 
divide by adding a suitable (compensating) fluid in the future cosmological scenario described 
by the models f l3.24p and f l3.25p . Here, we indicate a possible realization of it. We examine 
an inhomogeneous fluid with its energy density p, pressure P and the Eos parameter w as a 
function of p, i.e., u = uj{p). The EoS is expressed as 



z = -1. 




-0.85 



(a) 



(b) 




(4.17) 
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We explore the simple case 

a;(p) = Aoa(z)p'^^i-l, (4.18) 

where a is a constant and Aq is a positive parameter. Moreover, a{z) = —1 when z > 
and a{z) = 1 when z < 0, such that the fluid is in the phantom region for z > and in the 
quintessence region for 2; < 0. The fluid energy density reads 

p = Po {Bo - cx{z) \og{z + 1)) , (4.19) 

where po = [3(a — 1)Aq]^^^^~°'^ and Bq are positive parameters depending on the initial 
conditions. 

We note that one can choose Bq = 1 without the loss of generality and in this way 
the energy density is deflned as a positive quantity. If we take a > 1, when z — +00 
or 2 — > — l"*", the energy density asymptotically tends to zero. For z = 0, we have a 
maximum, p{z = 0) = Pq, so that we should require po ^ ^/'^^j namely, the fluid energy 
density is always small with respect to the dark energy density given by our models for the 
cosmological constant. A fluid in the form of Eq. fl4.19p may asymptotically produce a (Big 



to — tY , where t < to and (3 > 1 



for general study of singularities 



Rip) singularity H{t) ~ 
in modifled gravity, see j46|), only for (3 = l/{2a — 1) [47|, but in our case a > 1, so that 
this kind of divergence can never appear. If we add this fluid in the scenario described by 
F{R) gravity models in Eqs. (I3.24p and f l3.25p . when z — )■ — 1 we flnd 

^cfr = °\ ^ -1 + -f^ . 4.20 
Pde + P A/k^ 

This means that owing to the presence of fluid, the oscillations of the effective EoS parameter 
realize not around the phantom divide but around Ucs given by the last equation, namely 
in the quintessence region. With an accurate fltting of the parameters, in this way we may 
avoid the crossing of the phantom divide. 

We can also add a fluid to the cosmological scenario in order to have an asymptotical 
phantom phase without the Big Rip singularity. To this purpose, we investigate the EoS 
parameter of the fluid as in Eq. fl4.18p with Aq > and a{z) = —1, which describes a 
phantom fluid. The fluid energy density is given by Eq. f l4.19p . We put Bq = and a < 1 
such that 1/(1 — a) can be an even number and one can have the energy density deflned as 
a positive quantity. In this way, the fluid energy density decreases until 2 = and then it 
starts to grow up. We can take po sufficiently small so that the fluid contribution can become 
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dominant only in the asymptotical limit, when z is close to —1, avoiding the quintessence 
region in the final cosmological evolution of our F{K) gravity models. From the equation of 
motion ?>H'^ / h? = p, we obtain 



3 dz' , , 

In our case, it is easy to verify that t ~ | log(2 + l)(2a!-i)/(2a-2) | g^j^^ if a < 1/2, when 
z{t) — !■ —1 the integral diverges and t — )■ +oo, avoiding the Big Rip at a finite time. In this 



kind of models, the fiuid energy density increases with time^ 



Dut w — )■ — 1 asymptotically, so 



that there can be no future singularity. However, in Ref. |43| a careful investigation on the 
conditions necessary to produce this evolution has been done, and it has been demonstrated 
that this fiuid can rapidly expand in the future, leading to the disintegration of all bound 
structures (this is the so-called "Little Rip"). For example, a planet in an orbit of radius 
R around a star of mass M will become unbound when -(47r/3)(p + 3P)R^ ~ M. In 
our case, — (p + 3P) = Aqp"' and in the future every gravitationally bound system will be 
disintegrated |42|. 



C. Growth of the matter density perturbations: growth index 



In this subsection, we study the matter density perturbations. The equation that governs 
the evolution of the matter density perturbations for F{R) gravity has been derived in the 
literature (see, for example, 48| and references therein). Under the subhorizon approxima- 
tion, the matter density perturbation 6 = — satisfies the following equation: 



5 + 2H5 - 47rGefr(a, k)pra5 = 



(4.22) 



with k being the comoving wavenumber and Gcff(«, k) being the effective gravitational "con- 
stant" given by 



Gc{i{a,k) 



G 



F'{R) 



(P/a^) {F"iR)/F'iR)) 



(4.23) 



l + 3{kya^){F"{R)/F'{R))_ 
It is worth noting that the appearance of the comoving wavenumber k in the effective 
gravitational constant makes the evolution of the matter density perturbations dependent 
on the comoving wavenumber k. It can be checked easily, by taking F{R) = Rin Eq. (14.231) . 
that the evolution of the matter density perturbation does not have this kind of dependence 
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in the case of GR. In Fig. [TJ we show the cosmological evolution as a function of the redshift 
z and the scale dependence on the comoving wavenumber k of this effective gravitational 
constant for the case of model Fi{R) in Eq. fl3.24p . while in Fig. |H]we depict those for the case 
of model F2{R) in Eq. (13.251) . In both these cases, we have fixed 6 = 1 and used 7 = 1/1000. 
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FIG. 7: (a) Cosmological evolution as a function of z and the scale dependence on k of the effective 
gravitational constant G^s for the model Fi{R) with 6=1 and 7 = 1/1000. (b) Cosmological 
evolution of G^s as a function of z in the model Fi{R) with 6 = 1 and 7 = 1/1000 for k = IMpc^^ 
(blue), k = O.lMpc"^ (green), k = O.OlMpc"^ (red) and k = O.OOlMpc-^ (fuchsia). 
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FIG. 8: (a) Cosmological evolution as a function of z and the scale dependence on k of the effective 
gravitational constant G^s for the model F2{R) with 6=1 and 7 = 1/1000. (b) Cosmological 
evolution of G^q as a function of z for the model F2{R) with 6 = 1 and 7 = 1/1000. Legend is the 
same as Fig. [71 



Another important remark is to state that in deriving Eq. (I4.22p . we have assumed the 
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subhorizon approximation (see j49|). Namely, co moving wavelengths \ = a/k are considered 
to be much shorter than the Hubble radius as 

K^^H\ (4.24) 

This means that we examine the scales of \ogk > —3. On the other hand, as it was pointed 



out in Ref. |50[, for large k we have to take into account deviations from the linear regime. 
Hence, we do not consider the scales of log/c > — 1 and take the results obtained for logfc 
close to —1. 

Instead of solving Eq. f l4.22p for the matter density perturbation 6, we now introduce the 
growth rate /g = din 6 /din a and solve the equivalent equation to Eq. fl4.22p for the growth 
rate in terms of the redshift z, given by 

df.iz) ^ ri+^dHiz) _ 2 _ IM + 3 m^(l + ^)^Geff(a(z),fc) ^ ^ _ ^^^5) 



dz \H{z) dz •"'^'Jl + z 2 H^{z) G 

Unfortunately, Eq. (I4.25P cannot be solved analytically for the models Fi{R) and F2{R), 
but it can be solved numerically by imposing the initial conditions. Therefore, we execute 
the numerical calculations for both the model Fi{R) and the model F2{R) with the condition 
that at a very high redshift the growth rate becomes that in the ACDM model. In Fig. [HI we 
illustrate the cosmological evolution as a function of the redshift z and the scale dependence 
on the comoving wavenumber k of the growth rate for the model Fi{R), while we depict 
those of the growth rate for the model F2{R) in Fig. [TOl 

One way of characterizing the growth of the matter density perturbations could be to 
use the so-called growth index 7, which is defined as the quantity satisfying the following 
equation: 

f,{z) = fi^(^)^(^) , (4.26) 

with f2m(^) = ^3^^" being the matter density parameter. 

It is known that the growth index 7 in Eq. fl4.26p cannot be observed directly, but it can 
be determined from the observational data of both the growth factor fg{z) and the matter 
density parameter Qrai^) at the same redshift z. Even if the growth index is not directly 
observable quantity, it could have a fundamental importance in discriminating among the 
different cosmological models. One of the reasons is that in general, the growth factor 



fg{z), which can 



at different z 
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3e estimated from redshift space distortions in the galaxy power spectra 



52| . may not be expressed in terms of elementary functions and this fact 
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FIG. 9: (a) Cosmological evolution as a function of the redshift z and the scale dependence on the 
comoving wavenumber k of the growth rate /g for the model F\{K). (b) Cosmological evolution of 
the growth rate /g as a function of z in the model Fx{R) for k = O.lMpc"^ (green), k = O.OlMpc"^ 
(red) and k = O.OOlMpc"^ (blue). 
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FIG. 10: (a) Cosmological evolution as a function of the redshift z and the scale dependence on the 
comoving wavenumber k of the growth rate /g for the model F^iK). (b) Cosmological evolution of 
the growth rate /g as a function of z for the model F2 (R) ■ Legend is the same as Fig. [9l 



makes the comparison among the different models difficult. If Eq. fl4.26p is satisfied with 
any ansatz for the growth index 7, then its determination could provide an easy and fast 
way to distinguish between cosmological models. 

Various parameterizations for the growth index 7 have been proposed in the literature. In 
the first stage works on this topic, 7 was taken constant (see 53|). In the case of dark fiuids 
with the constant EoS uq in GR, it is 7 = 3 {uq — 1) / (6u;o — 5) (for the ACDM model, the 
growth index is 7 ~ 0.545). Although taking 7 constant is very appropriated for a wide class 



29 



ro 

0.45 
0.40 
0.35 
0.30 
0.25 
0.20 



ro 

0.50 ■ 

0.45 

0.40 

0.35 

0.30 

0.25 

0.20 



-2.5 



Logk(Mpc-') 



-2.5 -2.0 -1.5 -1.0 



,Logk(Mpc-') 



(a) (b) 

FIG. 11: Constant growth index as a function of log A; for the model Fi{R) (a) and for the model 
F2{R) (b). The bars express the 68% CL. 

of dark energy models in the framework of GR (for which |7'(0) | < 0.02), for modified gravity 
theories 7 is not constant in g eneral (the cases of some viable F{R) gravity models have been 
investigated in Refs. 50|, l5J]) and the measurement of |7'(0)| could be very important in 
order to discriminate between different theories. For this reason, another parameterizations 
has been proposed. The case of a linear dependence ■j{z) = ■Jq + 'JqZ was treated in Ref. [55 1. 
Recently, an ansatz of the type 7(2) =70 + 71^/(1 + 2;) with 70 and 71 being constants was 
explored in Ref. 56|and a generalization given by 7(z) = 70 + 7i;z/(l + z)'^ with a being 



a constant in Ref. 



50| . In the following, we study some of these parameterizations of the 



growth index for the case of the models Fi{R) and F2{R). 



1. 7 = 70 

We consider the ansatz for the growth index given by 

7 = 7o , 

where 70 is a constant. 

In Fig. [m we display the results obtained by fitting Eq. f l4.26p to the solution of Eq. fl4.25p 
for different values of the comoving wavenumber k for the two models Fi{R) and F2{R). We 
note that in these and following plots, the bars express the 68% confidence level (CL) and 
the point denotes the median value. The first important result for both models is that the 
value of the growth index has a strong dependence with log k. This scale dependence seems 
to be quite similar in both models. 
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In order to check the goodness of our fits, in Fig. [12] we show cosmological evolutions of 
the growth rate fg{z) and as functions of the redshift z together for several values 

of the comoving wavenumber k for the models Fi{R) and F2{R). To clarify these results, in 
Fig. [12] we also illustrate the cosmological evolution of the relative difference between fg{z) 
and Q^{zy° as a function of z for the same values of k in these models. The first remarkable 
thing is that for both models the function Q^{zy'^ fits the growth rate for large scales (i.e., 
lower k) very well, but this is not anymore the case for larger values of k. In fact, if we 
do not consider lower values for z (i.e., z < 0.2), for logk = —2 the relative difference is 
smaller than 3% for both models, while for log A; = — 1 can arrive up to almost 13%. For 
log A; = —3, we see that the relative difference is always smaller than 1.5% for the model 
Fi{R) and smaller than 1% for the model F2{R). 
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FIG. 12: Cosmological evolutions of the growth rate /g (red) and f^m (blue) with 7 = 70 as 
functions of the redshift z in the model Fi{R) for k = O.lMpc"^ (a), A; = O.OlMpc"^ (b) and 
k = O.OOlMpc^^ (c), and those in the model F2{R) for k = O.lMpc^^ (d), k = O.OlMpc"^ (e) and 
k = O.OOlMpc^^ (f). 
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FIG. 13: Cosmological evolution of the relative difference 



/g 



with 7 = 7o for k = O.lMpc 



(red), k = O.OlMpc ^ (blue) and k = O.OOlMpc ^ (green) in the model Fi{R) (a) and the model 
F2{R) (b). 

2- 7 = 70 + liz 

With the same procedure used in the previous subsection, we explore a hnear dependence 
for the growth index 



where 71 is a constant. 

In Fig. HU we depict the parameters 70 and 71 for several values of log A; in both the 
models. As is the same as the case 7 = 70, it can easily be seen that the scale dependence 
of the parameters 70 and 71 is similar in these models. We can also find that 70 ~ 0.46 
for the model Fi{R) when log/c < —2, whereas 70 ~ 0.51 for the model F2{R) when 
logk < —2.5. For both these models, the value of 71 has a strong dependence on k in the 
range of log A; > —2.25, but in the range of log /c < —2.25 this dependence becomes weaker. 

In Fig. [151 we illustrate cosmological evolutions of the growth rate fg{z) and flj^{zy^^'' 
as functions of the redshift z together for the models Fi{R) and F2{R). We can see that the 
fits for log A; = 0.1 have been improved in comparison with the same fits as the case with a 
constant growth index. Also, for logk < 0.1 the fits continue to be quite good. In order to 
demonstrate these facts quantitatively, in Fig. [16] we plot the cosmological evolution of the 
relative difference between fg{z) and ^ly^iz)"'^^^ as a function of z for several values of k in 
the models Fi{R) and F2{R). In this case, for log/c = —1 the relative difference is smaller 



7 



7o + 7iz , 



(4.27) 
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FIG. 14: Growth index fitting parameters in the case 7 = 70 + 7i-2 as a function of log A; for the 
model Fi{R) [(a) and (b)] and the model F2{R) [(c) and (d)]. Legend is the same as Fig. [TTl 

than 7.5% in both the models if we do not consider lower values for z (i.e., z < 0.2). We 
also see that the linear growth index improves the fits in both the models for log A; = —2 
in comparison with those for a constant growth index. In this case, the relative difference 
for the model Fi{R) is always smaller than 1%, whereas that for model F2{R) is smaller 
than 2%. Finally, for logk = —3 the results obtained for a constant growth index are quite 
similar to those for a linear dependence on z. 



Next, we examine the following ansatz for the growth index: 

7 = 70 + 717^- (4.28) 

1 + z 

In Fig. [T71 we depict the parameters 70 and 71 for several values of log k for both the 
models. The scale dependence of these parameters on k is shown. The behavior of the 
parameter 71 seems to be quite similar to that for the previous case 7 = 70 + 712;, but it is 
worth cautioning that the scale of the figures are different from each other, and that for the 
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(d) (e) (f) 

FIG. 15: Cosmological evolutions of the growth rate /g (red) and Q]n (blue) with 7 = 70 + 71^ 
as functions of the redshift z in the model Fi{R) for k = O-lMpc""*^ (a), k = O.OlMpc"^ (b) and 
k = O.OOlMpc^^ (c), and those in the model F2{R) for k = O.lMpc^^ (d), k = O.OlMpc"^ (e) and 
k = O.OOlMpc^^ (f). 
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FIG. 16: Cosmological evolution of the relative difference ' ^ ^ with 7 = 70 + 'Jiz for k = 
O.lMpc"^ (red), k = O.OlMpc"^ (blue) and k = O.OOlMpc^^ (green) in the model Fi{R) (a) and 
the model F2{R) (b). 

present ansatz the scale dependence of 71 is stronger than that for the previous case. It can 
also be seen that 70 ~ 0.465 for the model Fi{R) and 70 ~ 0.513 for the model F2{R) in the 
scale log A; < —2.5. 
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FIG. 17: 

the model 



Growth index fitting parameters in the case 7 = 70 + 7i as a function of log k for 
Fi{R) [(a) and (b)] and the model F2{R) [(c) and (d)]. Legend is the same as Fig. [TTJ 



In Fig. [TSl we plot cosmological evolutions of the growth rate fg{z) and VL^ (z)^(^) in the 
models Fi{R) and F2{R) for several values of k, as demonstrated in the previous subsections. 
We can see the fits for log k < —2 are quite good, as those in the previous ansatz for the 
growth index. In the case of higher values of log k, it seems that the fits are similar to those 
for a constant growth rate and these fits do not reach the goodness of those for the case of 
7 = 7o + liz- 

In order to analyze the fits quantitatively, in Fig. [T^ we display the cosmological evolution 
of the relative difference between fg{z) and Q^{zy^^^ for several values of k in the models 
Fi{R) and F2{R). We see that the relative difference for log/c = —1 is smaller than 12% 
(if we do not consider z < 0.2) for both the models. Thus, it is confirmed that these 
fits are better than those for the constant growth rate, but these are worse than those for 
7 = 7o + 7i^- For lower values of log k, the relative difference is smaller than 2% in 2; > 0.2. 

As a consequence, through the investigations of these different ansatz for the growth 
index, it is concluded that 7 = 70 + 71-2 is the parameterization that can fit Eq. fl4.26p to the 
solution of Eq. (14.251) better in a wide range of values for k. Even though the behavior of the 
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FIG. 18: Cosmological evolutions of the growth rate /g (red) and Q]a (blue) with 7 = 70 + lij^ 
as functions of the redshift z in the model Fi{R) for k = O.lMpc"^ (a), k = O.OlMpc"^ (b) and 
k = O.OOlMpc"^ (c), and those in the model F2{R) for k = O.lMpc"^ {d), k = O.OlMpc"^ (e) and 
k = O.OOlMpc^^ (f). 




(a) (b) 



FIG. 19: Cosmological evolution of the relative difference ' j ™' with 7 = 70 + lij^ for k = 
O.lMpc"^ (red), k = O.OlMpc"^ (blue) and k = O.OOlMpc"^ (green) in the model Fi{R) (a) and 
the model F2{R) (b). 

parameters 70 and 71 in the models Fi{R) and F2{R) is quite similar to each other, in order 
to distinguish between these models in Fig. [H] we can see that the more differences between 
these models come from the values of 70 for \ogk < —2. In fact, as remarked before, for 
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\ogk < —2.5 we have 70 ~ 0.46 for the model Fi{R) and 70 ~ 0.51 for the model F2{R). 



V. UNIFIED MODELS FOR EARLY AND LATE-TIME COSMIC ACCELERA- 
TION 

One of the most important goals for modified gravity theories is that all the processes 
in the expansion history of the universe are realized, namely, the universe starts with an 
inflationary epoch, followed by the radiation dominated era and the matter dominated uni- 
verse, and finally the late cosmic acceleration epoch is actually achieved without invoking 



the presence of dark components in the universe 
tion with dark energy, see Ref. 



11| (for a first F{R) theory unifying infla- 



17l|). Models of the type (13.21) may be combined in a natural 
way to obtain the phenomenological description of the infiationary epoch. For example, a 
'two-steps' model may be the smooth version, given by 

F{R) =R-2A[1- e-^/('^)] - Ai e{R - R-,) . (5.1) 

Here, 6{R — Rq) is the Heaviside's step distribution, Ri is the transition scalar curvature 
at inflationary scale and Ai is a suitable cosmological constant producing inflation, when 
R^ Ri. The main problem associated with this sharp model is the appearance of a possible 
antigravity regime in a region around the transition point between inflation and the universe 
described by the ACDM model. The antigravity in a past epoch is not phenomenologically 
acceptable. Furthermore, adding some terms would be necessary in order for inflation to 
end. 

In this section, we study two applications of exponential gravity to achieve an unifled 
description of the early-time inflation and the late-time cosmic acceleration. In particular, 
we show how it is possible to obtain inflationary universes with different numbers of e-folds 
by choosing different models parameters in the presence of ultrarelativistic matter in the 
early universe. 



Following the flrst proposal of Ref. [3J], we start with the form of F{R) with a natural 
possibility of a unifled description of our universe 

F{R) = R-2A{l-e-^)-A; 1 - e^^T + 7 (^-J-^^ i?" , (5.2) 

where Ri and Ai are the typical values of transition curvature and expected cosmological 
constant during inflation, respectively, and n is a natural number larger than unity (here, 
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we do not write the correction term for the stabihty of oscillations in the matter dominated 
era). In Eq. (15. 2p . the last term 7(l/i?"~^)i?°, where 7 is a positive dimensional constant 
and a is a real number, works at the inflation scale Ri and is actually necessary in order to 
realize an exit from inflation. 

We also propose another nice inflation model based on the good behavior of exponential 
function described as 



sm I Tre 

R 



F{R) = R-2A{l-e^m)-A, + ^ . (5.3) 



Here, the parameters have the same roles of the corresponding ones in the model in Eq. ( 15.20 . 
We note that the second term of the model vanishes when i? <^ i?j and tends to Ai when 
R ^ -Ri. We analyze these models, i.e.. Model I in Eq. (15. 2 p and Model II in Eq. (15. 3p . and 
explore the possibilities to reproduce the phenomenologically acceptable inflation. 



A. Inflation in exponential model (Model I) 



First, we investigate the model in Eq. (15.20 . For simplicity, we describe a part of it as 

1 



MR) = -AAl-e 



R". 



(5.4) 



We note that if > 1 and a > 1, when ^ i?i(~ RA, we obtain 



R"" 

h ^ 



(5.5) 



and the absence of the effects of inflation during the matter dominated era. We also find 



m) 



k;n{n - l)i?'^-2 
R^ 



+ 7« 



R 



a-l 



Ri 

^' +A^ 



(5.6) 



R? 



+ 7«(«-l)5^- (5.7) 



R? 



Since when R = Ri [(n — l)/n]^^" the negative term of f-'{R) has a minimum, in order to 
avoid the anti-gravity effects (this means, |//(-R)| < 1), it is sufficient to require 



,' n — 1\ " _ n-l 

Ri > AiU ( 1 e" 



(5.8) 



38 



It is necessary for the modification of gravity describing inflation not to have any influence 
on the stabihty of the matter dominated era in the small curvature limit. When R <^ R^, 
the second derivative of fl'{R), given by 



fRY-' _ , , fR^ 



-n(n-l)^-J +7a(a-l)^^ 



(5.9) 



must be positive, that is, 

n>a. (5.10) 

We require the existence of the de Sitter critical point -Rds which describes inflation in the 
high-curvature regime of /i(-R), so that fi{Rds ^ Ri) — —^i + tIV-^i""^) R"- this case, 
if we put -Ri = i?dS; we may solve Eq. fl2.4p for a constant curvature, and therefore we 
obtain (in vacuum, namely, if the effective modified gravity energy density is dominant over 
matter) , 

R.S = , . , f^V»l. (5.11) 



7(2 - a) + 1 ' V 

The last two conditions have to be satisfied simultaneously. By using Eq. (15. Sp . we also 
acquire 

fn-l\ " n-l 

>n{ e"— . (5.12) 



7(2 - a) + 1 \ n 
Instability and number of e- folds during inflation 

The well-known condition to have an instable de Sitter solution (see Sec. IIVBD is given 



by 

F'{Rds) ^ . TQ\ 

which leads to 

a7(a-2)>l, (5.14) 

for our model. Here, we have considered fi{Rds) — + 7(1/-^!°"^) R"- From Eqs. fl5.12p - 
(I5.14p . we have to require 

2 + l/7>a>2. (5.15) 
Thus, we may evaluate the characteristic number of e-folds during inflation 

iV = log^, (5.16) 

Ze + 1 
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where Zi and are the redshifts at the beginning and at the end of early time cosmic 
acceleration. Given a small cosmological perturbation yi{z]) at the redshift 2i, we have from 
Eq. fl4.13p avoiding the matter contribution 



y,{z,) = C,{z, + l)\ 



(5.17) 



with 



X = 




25 



) 



(5.18) 



where a; < if the de Sitter point is unstable. Thus, the perturbation yi{z) in Eq. fl4.13p 
grows up in expanding universe as 



Here, we have considered Cq = yi{z{)/{zi + 1)^. When yi{z) is on the same order of the 
effective modified gravity energy density y^ of the de Sitter solution describing inflation (we 
remind, yo = -Rds/(12?7i^)), the model exits from inflation. We can estimate the number of 
e-folds during inflation as 



A value demanded in most inflationary scenarios is at least = 50-60. 

A classical perturbation on the (vacuum) de Sitter solution may be given by the presence 
of ultrarelativistic matter in the early universe. The system gives rise to the de Sitter solution 
where the universe expands in an accelerating way but, suddenly, it exits from inflation and 
tends towards the minimal attractor at i? = (the trivial de Sitter point). In this way, the 
small curvature regime arises and the physics of the ACDM model is reproduced. 

B. Inflation in Model II 

Next, we study the inflation model in Eq. fl5.3p . By performing a similar analysis to that 
in the previous subsection, we find that also in this case, if a > 1 and n > 1, we avoid the 
effects of inflation at small curvatures and it does not influence the stability of the matter 
dominated era. The de Sitter point exists if = R^s and it reads as in Eq. (15. lip under 
the condition [R/Ri)"' ^ 1. Thus, the inflation is unstable if the condition in Eq. f l5.14p 
is satisfied. The bigger difference between the two models exists in those behaviors in the 



yi{z) = yi{zi) 



(5.19) 




(5.20) 
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transition phase between the small curvature region (where the physics of the ACDM model 
emerges) and the high curvature region. This means that the no antigravity condition is 
different in the two models and such a condition becomes more critical in the transition 
region. Therefore, in the following we are able to make the different choices of parameters in 
the two models. We note that since dark energy sector of the above models only originates 
from exponential gravity, all qualitative results in terms of the behavior of the dark energy 
component in exponential gravity found in the previous sections remain to be valid. 



VI. ANALYSIS OF INFLATION 



In this section, we perform the numerical analysis of the early time acceleration for the 
unified models in Eqs. (15. 2p and (15. 3p . by choosing appropriate parameters according with 
the analysis in Sec. |Vl For this aim, it is worth rewriting Eq. ( 12.13^ by introducing a suitable 
scale factor at the inflation. We can choose = Ai. The effective modified gravity 
energy density yniz) is now defined as 

Here, we have neglected the contribution of standard matter and supposed the presence of 
ultrarelativistic matter/radiation in the hot universe scenario, whose energy density prad at 
the redshift equal to zero is related with the scale as 

X = ^ • (6-2) 

Since the results are independent of the redshift scale, we set 2; = at some times around 
the end of inflation. Equation (I2.13p reads 



+ 



z+l\ 2M^F"{R) [yniz] + x{z + 1)^] 

Mz) 2 - F'{R) 

[z + 1)2 M^F"{R) [yH{z) + x{z + l\ 
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, {F'{R) - l)2x{z + 1)^ + {F{R) - R)/M^ ^ 
{z + lY2M^F"{R)[yH{z)+x{z + lY] ' 



Moreover, the Ricci scalar is expressed as 

dyu^z) 



R = M' 



Ayniz) -{z + 1) 



dz 



(6.4) 
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Thus, it is easy to verify tliat in tfie de Sitter universe witli R = Rds the perturbation yi{z) on 
the solution z/q = i?ds/(4 M^) is effectively given by Eq. (15.171) . i.e., yi = Co{z+iy, according 
with Eq. (14.131) if we neglect the contribution of standard matter. In this derivation, we have 
assumed the contribute of ultrarelativistic matter to be much smaller than y^. However, as 
stated above, this small energy contribution may originate from the perturbation yi{zi) at 
the beginning of inflation, which, if x < 0, grows up in the expanding universe making 
inflation unstable. 



Model I 

First, we explore the model in Eq. (15. 2p . We have to choose the parameters as Ai ~ 
]^Qioo-i20yY^ The dynamics of the system is independent of this choice. Here, we summarize 
the conditions for inflation already stated in Sec. IV AI 

n-l 

\ 71 71 — 1 

j e~^^ , (no antigravity effects) 

-Ri = -Rds 5 07(0; — 2) > 1 , ( ) ^ ^ ' (existence of unstable dS solution) 

with 77, > 1, 2 + 1/7 > a > 2 and -Rds = 2Ai/ [7(2 — a) + 1]. Since 7 and a are combined in 
7(a — 2), we can fix 7 = 1, so that -Rds = 2Ai/ (3 — a) and 3 > a > 2. The instability factor 
X in Eq. (I5.18P only depends on i?ds- Hence, by studying the phenomenology of inflation, 
we examine the variation of a parameter (and, as a consequence, that of Ri). We take 
n = 4 and Ri = 2Ai, which satisfy the condition for no antigravity well. We analyze three 
different cases of a = 5/2, 8/3, and 11/4. In these cases, we have -Rds = 4Ai, 6Ai, and SAi, 
respectively. 

In Fig. [201 we plot the cosmological evolution of the quantity 2F{R/Ai) - {R/ K-^F' {R/ K,) 
as a function of the redshift z in the three cases. The value of "zero" of this quantity 
corresponds to the de Sitter points of the model. We can recognize the unstable de Sitter 
solutions of inflation and the attractor in zero (the de Sitter point of current acceleration is 
out of scale). 

Despite the fact that the three considered values of a are very close each other, the values 
of i?ds and X significantly change and the reactions of the system to small perturbations are 
completely different. By starting from Eq. (I5.20p . we may reconstruct the rate yi{zi)/yo 
between the abundances of ultrarelativistic matter /radiation and modified gravity energy at 
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(a) (b) (c) 

FIG. 20: Cosmological evolution of the quantity 2F{R/h.\) — {R/Ai)F'{R/Ai) as a function of the 
redshift z for exponential model with a = 5/2 (a), a = 8/3 (b) and a = 11/4 (c). The "zeros" of 
these graphics indicate the de Sitter solutions of the model. 

the beginning of inflation in order to obtain a determined number of e-folds during inflation 
in the three different cases, by taking into account that x = —0.086, —0.218, and —0.270 for 
a = 5/2, 8/3, and 11/4, respectively. For example, in order to have = 70, for a = 5/2, a 
perturbation of yi(2;i)/?/o ~ 10~^ is necessary; for a = 8/3, a perturbation of ?/i(2;i)/i/o ~ 10"''' 
is sufficient; whereas for a = 11/4, yi{zi)/yo ~ 10~^. The system becomes more unstable, 
as (3 — a) is closer to zero. 

In studying the behavior of the cosmic evolution in Model I for the three different cases, 
we set X = 10~^yo/(^i + 1)^ in Eq. (16. 3 p for the case a = 5/2 and x = ^^"^Uo/i^i + 
1)^ for the cases a = 8/3, 11/4. In these choices, the effective energy density originating 
from the modification of gravity is 10^ and 10^ times larger than that of ultrarelativistic 
matter/radiation during inflation. By using Eq. fl5.20p . we can predict the following numbers 
of e-folds: 

N ~ 107 (for a = 5/2) , 
iV ~ 64 (for a = 8/3) , 

^ 51 (for a = 11/4) . (6.5) 

In order to solve Eq. (16. 3 p numerically, we use the initial conditions 

dyniz) 



d{z) 



yH[z) 







(6.6) 



4Ai ' 

at the redshift 2i ^ when inflation starts. We put zi = lO'^^, 10^'', and 10^^ for a = 5/2, 
8/3, and 11/4, respectively (just for a more comfortable reading of the graphics). We also 
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remark that the initial conditions are subject to an artificial error that we can estimate to 
be in the order of exp [— (i?ds/-Ri)"] ~ 10^''. This is the reason for which we only consider 




(d) (e) (f) 

FIG. 21: Plots of yn [a-c] and wmg [d-f] as functions of the redshift z for Model I with a = 5/2 
[a-d], a = 8/3 [b-e] and a = 11/4 [c-f]. 

In Fig. [211 we illustrate the cosmological evolutions of yn and the corresponding modified 
gravity EoS parameter umg (defined as in Eq. f l3.15p ) as functions of the redshift z in the 
three cases. We can see, during infiation umg is indistinguishable from the value of -1 and 
Hh tends to decrease very slowly with respect to yn = 1,3/2,2 for a = 5/2,8/3, 11/4, so 
that the curvature can be the expected de Sitter one, Rdsi= ^Uh) = 4Ai,6Ai,8Ai. The 
expected values of at the end of infiation may be derived from the number of e-folds in 
(16. 5p during infiation and read ~ —0.47 for a = 5/2; z^ — —0.74 for a = 8/3; Zc — —0.39 
for a = 11/4. The numerical extrapolation yields 



Vh 


>e) 


= 0.83yH 






= 0.825i2ds , 


(for 


a 


= 5/2) 


Vh 


>e) 


= 0.88yH 


>i) , 


R{Ze) 


= 0.853i?ds , 


(for 


a 


= 8/3) 


Vh 


>e) 


= 0.92yH 


>i) , 


R{Ze) 


= 0.911/2ds. 


(for 


a 


= 11/4) 



To confirm the exit from infiation, in Fig. [22] we plot the cosmological evolutions of yn and 
R/Ai as functions of the redshift z in the region —l<z<l, where Ze is included. The 
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effective modified gravity energy density and the curvature decrease at the end of inflation 
and the physical processes described by the ACDM model can appear. 




Model II 

Next, we investigate Model II in Eq. (15. 3p . Here, in order to satisfy the condition for 
no antigravity we choose n = 3 and Ri = 2A[, so that F'{R > 0) > 0. We take 7 = 1 
again and we execute the same numerical evaluation for a = 5/2, 13/5,21/8 in this model 
as that in the previous case for Model I. The corresponding de Sitter curvatures of inflation 
are i?ds = 4Ai,5Ai, 16Ai/3. Now, we obtain the factor in Eq. fl5.18p for instability as x = 
—0.086, —0.170, and —0.188 for a = 5/2, 13/5, and 21/8, respectively. Hence, we set 
X = 10'3^o/(^i+l)^fora = 5/2, x = 10"^ yo/l^i+l)^ for a = 13/5, and x = 10"^ yo/(^i+l)^ 
for a = 21/8. As a consequence, the numbers of e- folds during inflation result in = 80, 
54, and 61. The initial conditions are the same as those in the previous case in ( 16. 6p . 
Furthermore, we put Zi = W^^, 10^^, and 10^^ for a = 5/2, 13/5, and 21/8. 

Through the numerical extrapolation, we acquire the expected values of at the end of 
inflation as z^ = —0.80, —0.97, and —0.71, and the following values for the effective modified 
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gravity energy density and the Ricci scalar: 



ynize) = 0.82yH{z,) , R{z,) = 0.813i?ds , (for a = 5/2) 
yniz,) = OMynizi) , R{z,) = 0.884i?ds , (for a = 13/5) 
yniz,) = 0.79yH{z,) , R{z,) = 0.780/2ds • (for a = 21/8) 

For this model, in Fig. [23] we depict the cosmological evolutions of yn and R/Ai as functions 
of the redshift z in the region — 1 < 2; < 1 at the end of inflation. Again in this case, the 
effective modified gravity energy density and curvature decrease, and therefore inflation ends 
and then the physical processes described by the ACDM model can be realized. 




(d) (e) (f) 

FIG. 23: Cosmological evolution of yn [a-c] and R/Ai [d-f] as functions of the redshift z in the 
region -1 < z < 1 for Model II with a = 5/2 [a-d], a = 13/5 [b-e] and a = 21/8 [c-f]. 



VII. CONCLUSIONS AND GENERAL REMARKS 

In the present paper, we have examined a generic feature of viable F{R) gravity models, 
in particular, exponential gravity as well as a power form model. We have shown that the 
behavior of higher derivatives of the Hubble parameter may be affected by large frequency 
oscillations of effective dark energy, and consequently solutions may become singular and un- 
physical at a high redshift. The analyzed models approach to a model with the cosmological 
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constant in a manner different from each other, and hence it is reasonable to expect that the 
found results can be generalized to realistic F{R) gravity models, in which the cosmological 
evolutions are similar to those in a model with the cosmological constant. To support our 
claim, in the first part of this paper we have explicitly demonstrated how the origin of the 
problem influences the stability conditions satisfied by these models in order to reproduce 
the realistic matter dominated era. Since the corrections to the Einstein equations at the 
small curvature regime may lead to undesired effects at the high curvature regime, we have 
reconstructed a correcting (compensating) term added to the models in order to stabilize 
the oscillations of the effective dark energy in the matter dominated era with retaining the 
viability properties. It is emphasized that all the results we have found in an analytical way 
via studying the perturbation theory are confirmed by the numerical analysis performed on 
the models under consideration. Moreover, a detailed investigation on the cosmological evo- 
lutions of the universe described by those models has been executed. In particular, we have 
demonstrated that our correction term does not cause any problem to the viability of the 
models, and that the obtained results are consistent with recent very accurate observational 
data of our current universe and easily pass the local tests of the solar system. Furthermore, 
we have shown that the effective crossing of the phantom divide, which characterizes the de 
Sitter epoch, occurs in the very far future. A way to avoid the crossing of phantom divide 
by using inhomogeneous fluids has also been explored. 

After the discovery of the accelerated expansion of our universe, a lot of theories are 
proposed in order to explain it. The issue of discriminating among all of these theories has 
become very important. The first step in order to distinguish between theories can be the 
study of their expansion history, but it has been revealed that sometimes different models 
exhibit the same (or very similar) expansion history. For this reason, the investigation of 
growth of the matter density perturbations by using the so-called growth index can provide 
a significant tool in order to distinguish among the different gravitational theories. In this 
context, the growth of the matter density perturbations has been examined for our models. 
Several ansatz for the growth index have been considered, and consequently it has been 
concluded that the choice of the growth index as 7 = 70 + 712; is the most appropriate 
parameterization for these theories. 

In addition, in the second part of this paper we have discussed the inflationary cosmology 
in two exponential gravity models. It has explicitly been shown that different numbers of 
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e-folds during inflation can be obtained by taking different model parameters in the presence 
of ultrarelativistic matter, the existence of which makes inflation to end and realize the exit 
from it. We have performed the numerical analysis of the inflationary stage in two viable 
exponential gravity models. It has been found that at the end of the inflation, the effective 
energy density and therefore the curvature of the universe become small. As a result, we 
have proved that it is possible to acquire a gravitational alternative scenario for a unified 
description of inflation in the early universe with the late-time cosmic acceleration due to 
the ACDM-like dark energy domination. 

Finally, we mention that as an important future work in terms of our present investiga- 
tions, at the next step we plan to study cosmological perturbations 
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in such resultant 



F{R) gravity theories. We calculate the power spectrum of the cosmological perturbations 
as well as the tensor-to-scalar ratio in these models and compare those with the observational 
data from such as WMAP satellite future PLANCK satellite 59, 60|, QUIET 61, 62|, 



B-Pol 63[ and LiteBIRD 6J] in terms of the polarization of the CMB radiation. Further- 
more, it is meaningful to remark that the growth of the matter density perturbations in 
modified gravity affects the spectrum of weak lensing (for a concrete way of comparing the 



theoretical predictions with the observations, see [66|), and therefore more precise future 
observations of weak lensing effects have a potential to present the chance to find out the 
signal of the modification of gravity. 
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Appendix A: Conformal transformation of exponential model for inflation 

In several suitable conformal frame to study inflation may be the so-called "Ein- 

stein frame" . An F{R) gravity theory can be rewritten in the scalar field theory form via the 
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conformal transformation. We can rewrite the action in Eq. (12. ID by introducing a scalar 
field which couples to the curvature. Of course, this is not exactly physically-equivalent 
formulation, but the formulation in the Einstein frame may be used to obtain some of in- 
termediate results in simpler form (especially, the case that the matter is not taken into 
account). 

We introduce a scalar field A into the action 



/JF = ^ / [F'{A) {R-A) + F{A)] d'x . 



(7.1) 



Here, the subscript "JF" means "the Jordan frame" and we neglect the contribute of matter. 
By making the variation of the action with respect to A, we have A = R. We define the 
scalar field a as 



V2a 



\n[F'{A)] . 



We make the conformal transformation of the metric 



(7.2) 



for which we acquire the "Einstein frame" (EF) action of the scalar field a 65| 



'EF 



d xW—g 



R 1 [F"{A) 



M 



r'^dMA 



d —g 



M 



2^2 2 V^'(^) 
_^ - ^~g^-^d,ad^a + Via) 



A 



2k^ \F'{A) ' F'{Af 



where 



V{a) 



A 



F{A) 



(7.3) 



(7.4) 



{e"i?(e-") -e2"F[i?(e-")]} . (7.5) 



2^2 \F'{A) F'{Ay J 2k2 

Here, R{e~") is the solution of Eq. (17. 2p with A = R, becoming R a function of , and R 
denotes the Ricci scalar evaluated with respect to the conformal metric g^^. Furthermore, 
g = e~^°"(y' is the determinant of conformal metric. 

As an example, we explore our unified model (15. 2 p with 7 = 1. Since we are interested 
in the de Sitter solution, we take exp[—{R/Ri)"'] — )■ and neglect the cosmological constant 
A. In this case the potential V{a) reads 



Via) 



R 



a 



(e" - 26^") + A; e^^ 



(7.6) 
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According with Sec. IV A| we put Ri = Rds- It is clearly seen that for R = R^s, Cds 



— a/3/ (2^2) log(l + a) and V'{ads) = 0, where the prime denotes the derivative with respect 
to the inflation fleld a. Since ^"(0-^3) > 0, the scalar potential has a minimum, that is a 
necessary condition for a slow-roll inflation. For slow-roll parameters, we have to require 



1 

2^ 



V'(a) 



V"(a) 



< 1 



Via) 



< 1. 



(7.7) 



By deflning the energy density and pressure of cr as p^- = /2 — V{o-) and = a'^ /2 + V{a), 
these conditions imply that the gravitational fleld equations in the flat FLRW space-time 
are given by 3H^/k,^ = —V{a), 3H& ^ —V'{a), and that d{t) > 0, and hence guarantee a 
sufficiently long time inflation. In our case, since V^(o"ds) 7^ 0, these two conditions are well 
satisfled around the de Sitter solution. Thus, since d ~ 0, we find H^s = Rds/ [12(1 + a)] = 
^ds/12. 



Appendix B: Asymptotically phantom or quintessence modified gravity 

In general, realistic models of modifled gravity are similar to GR with the cosmolog- 
ical constant, i.e., the dark energy fluid with the EoS parameter cude = —1 and the 
de Sitter universe as the flnal scenario for the cosmological evolution. Since in principle 
quintessence/phantom-dark energy phases are not excluded by observations, it may be of 
some interest to try to reconstruct an F{R) gravity theory where the quintessence or phan- 
tom dark energy (with a constant wde) emerges. The big difficulty is due to the fact that in 
the dark energy density pde and pressure Pde of modifled gravity, the effective gravitational 
terms appear. In this appendix, we reconstruct the form of F{R) gravity which resembles 
to a fluid with wde being very close but not equal to —1. 

If the energy density of a quintessence / phantom fluid is given by 

p = poiz + (7.8) 
where cu is the EoS parameter, the Hubble parameter reads 

H{z) = ^ + \^/3^oil + uj) \og\z + 1] . (7.9) 
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Here, we have taken into account that u is very close to —1. We can write i? as a function 
of the redshift as 



= li^^Po [2 + 3(1 + co) \og{z + 1)] [1-3LU + 6(1 + cu) \og{z + 1)] 



(7.10) 



In addition, from Eq. fl2.10p . in vacuum we find 



Poff = Pde 



2^2 



dR{z)\^ dF{z) 



dz 

dR{z 
dz 



dz 



R{z) - F{z) 



6H^(z) 



dR{z)Y^ dF{z) 



dz 



dz 



d^R{z) \ dF{z) fdR{z) \ d^F{ 



dz^ 



dz 



dz 



dz^ 



(7.11) 



Here, F{R) model is expressed as a function of the redshift F{z). By equating p^s to p of 
Eq. (17.81) . we can find the F{R) model realizing this cosmology. For — 1| ^0, the solution 
of Eq. (17. lip is given by 

6kVo [11 + (34 - 9uj)lu] 



F z) 



From Eq. f l7.10p . we have 



(5 - 3ujy 



6k>o(1 + w) log(z + 1) 



(7.12) 



z = —1 + exp 



pon'^ib - 3uj) ± (1 + w) v/p;;7?]T6^T9A^^?(rTwH 

12poK^{l+Uj) 



(7.13) 



where the plus sign corresponds to the quintessence solution, whereas the minus sign does 
to the phantom one. We can now write the modified gravity model as a function of the Ricci 
scalar as 



n^ Po/^'(257 + 183u; + 27^2 _ 37,^3) ^ ^^^2 [I6i? + 9po«:'(l + w)'] 
r [R) = — — ^^1^ ± 



(7.14) 



2(5-3a;)2 2 

where po is a free parameter of the theory, and u is the EoS parameter of dark energy coming 
from the modification of gravity and equivalent to cjde- In this way, we have reconstructed 
the form of F{R) gravity that gives the quintessence or phantom fiuid solution in the empty 
universe. Remind that this reconstruction is valid for cjde close to —1. 
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